MV-based relative electron density estimation (iMREDe) for MR-LINAC dose calculation.

BACKGROUND: MR-LINAC systems have been increasingly utilized for real-time imaging in adaptive treatments worldwide. Challenges in MR representation of air cavities and subsequent estimation of electron density maps impede planning efficiency and may lead to dose calculation uncertainties. PURPOSE: To demonstrate the generation of accurate electron density maps using the primary MV beam with a flat-panel imager. METHODS: The ViewRay MRIdian MR-LINAC system was modeled digitally for Monte Carlo simulations. Iron shimming, the magnetic field, and the proposed flat panel detector were included in the model. The effect of the magnetic field on the detector response was investigated. Acquisition of projections over 360 degrees was simulated for digital phantoms of the Catphan 505 phantom and a patient treated for Head and Neck cancer. Shim patterns on the projections were removed and detector noise linearity was assessed. Electron density maps were generated for the digital patient phantom using the flat-panel detector and compared with actual treatment planning CT generated electron density maps of the same patient. RESULTS: The effect of the magnetic field on the detector point-spread function (PSF) was found to be substantial for field strengths above 50 mT. Shims correction in the projection images using air normalization and in-painting effectively removed reconstruction artifacts without affecting noise linearity. The relative difference between reconstructed electron density maps from the proposed method and electron density maps generated from the treatment planning CT was 11% on average along all slices included in the iMREDe reconstruction. CONCLUSIONS: The proposed iMREDe technique demonstrated the feasibility of generating accurate electron densities for the ViewRay MRIdian MR-LINAC system with a flat-panel imager and the primary MV beam. This work is a step towards reducing the time and effort required for adaptive radiotherapy in the current ViewRay MR-LINAC systems.

Monte Carlo model of a prototype flat-panel detector for multi-energy applications in radiotherapy.

BACKGROUND: The incorporation of multi-energy capabilities into radiotherapy flat-panel detectors offers advantages including enhanced soft tissue visualization by reduction of signal from overlapping anatomy such as bone in 2D image projections; creation of virtual monoenergetic images for 3D contrast enhancement, metal artefact reduction and direct acquisition of relative electron density. A novel dual-layer on-board imager offering dual energy processing capabilities is being designed. As opposed to other dual-energy implementation techniques which require separate acquisition with two different x-ray spectra, the dual-layer detector design enables simultaneous acquisition of high and low energy images with a single exposure. A computational framework is required to optimize the design parameters and evaluate detector performance for specific clinical applications. PURPOSE: In this study, we report on the development of a Monte Carlo (MC) model of the imager including model validation. METHODS: The stack-up of the dual-layer imager (DLI) was implemented in GEANT4 Application for Tomographic Emission (GATE). The DLI model has an active area of 43×43 cm2 , with top and bottom Cesium Iodide (CsI) scintillators of 600 and 800 µm thickness, respectively. Measurement of spatial resolution and imaging of dedicated multi-material dual-energy (DE) phantoms were used to validate the model. The modulation transfer function (MTF) of the detector was calculated for a 120 kVp x-ray spectrum using a 0.5 mm thick tantalum edge rotated by 2.5o . For imaging validation, the DE phantom was imaged using a 140 kVp x-ray spectrum. For both validation simulations, corresponding measurements were done using an initial prototype of the imager. Agreement between simulations and measurement was assessed using normalized root mean square error (NRMSE) and 1D profile difference for the MTF and phantom images respectively. Further comparison between measurement and simulation was made using virtual monoenergetic images (VMIs) generated from basis material images derived using precomputed look-up tables. RESULTS: The MTF of the bottom layer of the dual-layer model shows values decreasing more quickly with spatial frequency, compared to the top layer, due to the thicker bottom scintillator thickness and scatter from the top layer. A comparison with measurement shows NRMSE of 0.013 and 0.015 as well as identical MTF50 of 0.8 mm1 and 1.0 mm1 for the top and bottom layer respectively. For the DE imaging of the DE-phantom, although a maximum deviation of 3.3% is observed for the 10 mm aluminum and Teflon inserts at the top layer, the agreement for all other inserts is less than 2.2% of the measured value at both layers. Material decomposition of simulated scatter-free DE images gives an average accuracy in PMMA and aluminum composition of 4.9% and 10.3% for 11-30 mm PMMA and 1-10 mm aluminum objects respectively. A comparison of decomposed values using scatter containing measured and simulated DE images shows good agreement within statistical uncertainty. CONCLUSION: Validation using both MTF and phantom imaging shows good agreement between simulation and measurements. With the present configuration of the digital prototype, the model can generate material decomposed images and virtual monoenergetic images.

The effect of common dental fixtures on treatment planning and delivery for head and neck intensity modulated proton therapy.

PURPOSE: Proton treatment plan perturbation by common dental fixtures such as amalgams (Am) and porcelain-fused-to-metal (PFM) crowns has, to date, been uncharacterized. Previous studies have been conducted to determine the physical effect of these materials within the beam path for single spots, but their effects on complex treatment plans and clinical anatomy have not yet been quantified. The present manuscript aims to study the effect of Am and PFM fixtures on proton treatment planning in a clinical setting. METHODS: An anthropomorphic phantom with removable tongue, maxilla, and mandible modules was simulated on a clinical computed tomography (CT) scanner. Spare maxilla modules were modified to include either a 1.5 mm depth central groove occlusal amalgam (Am) or a porcelain-fused-to-metal (PFM) crown, implanted on the first right molar. Modified tongue modules were 3D printed to accommodate several axial or sagittal oriented pieces of EBT-3 film. Clinically representative spot-scanning proton plans were generated in Eclipse v.15.6 using the proton convolution superposition (PCS) algorithm v.15.6.06 using a multi-field optimization (MFO) technique with the goal of delivering a uniform 54 Gy dose to a clinical target volume (CTV) typical of a base-of-tongue (BoT) treatment. A typical geometric beam arrangement of two anterior oblique (AO) beams and a posterior beam was employed. Plans optimized without any material overrides were delivered to the phantom A) without implants; B) with Am fixture; or C) with PFM crown. Plans were also reoptimized and delivered with inclusion of material overrides to equate relative stopping power of the fixture with that of a previously measured result. RESULTS: Plans exhibit slightly greater dose weight towards AO beams. The optimizer accounted for inclusion of fixture overrides by increasing beam weights to the beam closest to the implant. Film measurements exhibited cold spots directly within the beam path through the fixture in plans with and without overridden materials. Cold spots were somewhat mitigated in plans including overridden materials in the structure set but were not entirely eliminated. Cold spots associated with Am and PFM fixtures were quantified at 17% and 14% for plans without overrides, respectively, and 11% and 9% with using Monte Carlo simulation. Compared with film measurements and Monte Carlo simulation, the treatment planning system underestimates the dose shadowing effect in plans including material overrides. CONCLUSIONS: Dental fixtures create a dose shadowing effect directly in line with the beam path through the material. This cold spot is partially mitigated by overriding the material to measured relative stopping powers. Due to uncertainties in modeling perturbation through the fixture, the magnitude of the cold spot is underestimated using the institutional TPS when compared to measurement and MC simulation.

Physical characterization of therapeutic proton delivery through common dental materials.

PURPOSE: Dental fixtures are commonplace in an aging, radiation treatment population. The current, local standard of practice in particle therapy is to employ treatment geometries to avoid delivery through implanted dental fixtures. The present study aims to observe the physical effect of delivering therapeutic proton beams through common dental fixture materials as prelude to an eventual goal of assessing the feasibility of using treatment geometries not specified for avoidance of oral implants. A sampling of common dental materials was selected based on prosthodontic consult and was evaluated in terms of relative stopping power and three-dimensional (3D) dose perturbation. METHODS: Amalgams, porcelain-fused-to-metal (PFM) crowns consisting of zirconia and non-noble base metals, and lithium disilicate implants were chosen for analysis. Theoretical stopping power (S) and mass stopping power (S/ρ) were calculated using the Stopping and Range of Ions in Matter (SRIM) application, basing stoichiometric compositions of each fixture on published materials data. S and S/ρ were calculated for a range of historically available compositions of amalgams from 1900 until the current era. The perturbance of S and S/ρ as a function of clinically relevant ranges of amalgam compositions for the modern era was analyzed. Water equivalent thickness (WET) and relative stopping power (Srel ) of each material was measured for a clinical spot-scanning proton beam with monoenergies of 159.9 and 228.8 MeV with a multi-layer ionization chamber (MLIC). Subsequently, 3D dose perturbation was assessed by delivering proton beams through a custom phantom designed to simulate both en-face and on-edge treatment geometries through the selected materials. A treatment plan mimicking the experimental delivery was constructed in the institutional treatment planning system and calculated using TOPAS-based Monte Carlo simulation (MCS). Experimental results were used to validate the MCS. Finally, treatment planning system (TPS) outputs were compared to MCS to determine the accuracy of the dose calculation model. RESULTS: Historical compositions of amalgams ranged in S from 44.8 to 42.9 MeV/cm, with the greatest deviation being observed for the 1900-1959 era. Deviation as a function of amalgam composition from the modern era was most sensitive to proportion of Hg, accounting for deviations up to -4.2% at the greatest clinically relevant concentration. S/ρ was not found to vary greatly between each porcelain and metal alloy material for PFM type crowns. Relative stopping powers ranged between 1.3 and 5.4 for all studied materials, suggesting substantial changes in proton range with respect to water. Film measurements of pristine spots confirm dose perturbance and shortening of proton range, with an upstream shift of each Bragg peak being observed directly behind the installed fixture. At high energies, cold spots were found in all cases directly behind each material feature with a medial fill-in of dose occurring distally. Qualitative agreement of spot perturbance was confirmed between film measurements and MCS. Finally, when comparing integrated depth doses (IDD) by summing over all axial directions, good agreement is observed between TPS and MCS. CONCLUSIONS: All dental materials studied substantially perturbed the dosimetry of pristine proton spots both in terms of WET/Srel as well as the spatial distribution of dose. Proton range was quantifiably shortened, and each dental material affected a cold spot directly behind the object with medial dose back-filling was observed distally. MCS and Eclipse dose calculations exhibited good agreement with measurements, suggesting that treatment planning without employing avoidance strategies may be possible with further investigation.

A novel method for fast image simulation of flat panel detectors.

We have developed a novel method for fast image simulation of flat panel detectors, based on the photon energy deposition efficiency and the optical spread function (OSF). The proposed method, FastEPID, determines the photon detection using photon energy deposition and replaces particle transport within the detector with precalculated OSFs. The FastEPID results are validated against experimental measurement and conventional Monte Carlo simulation in terms of modulation transfer function (MTF), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), contrast, and relative difference of pixel value, obtained with a slanted slit image, Las Vegas phantom, and anthropomorphic pelvis phantom. Excellent agreement is observed between simulation and measurement in all cases. Without degrading image quality, the FastEPID method is capable of reducing simulation time up to a factor of 150. Multiple applications, such as imager design optimization for planar and volumetric imaging, are expected to benefit from the implementation of the FastEPID method.

A Monte Carlo study of the impact of phosphor optical properties on EPID imaging performance.

We have developed a Monte Carlo computational model of a clinically employed electronic portal imaging device (EPID), and demonstrated the impact of phosphor optical properties on the imaging performance. The EPID model was built with Geant4 application for tomographic emission. Both radiative and optical transport were included in the model. Modulation transfer function (MTF), normalized noise-power spectrum times the incident x-ray fluence (qNNPS), and detective quantum efficiency (DQE) were calculated for simulated and measured data, and their agreement was quantified by the normalized root-mean-square error (NRMSE). MTF was computed using a 100 µm wide slit tilted by 1.5° and qNNPS was estimated using the Fujita-Lubberts-Swank method. DQE was calculated from MTF and qNNPS data. The NRMSE value was 0.0467 for MTF, 0.0217 for qNNPS, and 0.0885 for DQE, showing good agreement between measurement and simulation. Five major optical properties, phosphor grain size, phosphor thickness, phosphor refractive index, binder refractive index, and packing ratio were tested for their influence on the qNNPS, MTF, and DQE(0) of the model. Generally, the effect on the qNNPS is greater than MTF, and no impact on DQE(0), except from phosphor thickness, was observed. Multiple applications, such as imager design optimization and investigations of the dosimetric performance, are expected to benefit from the validated model.

Physics considerations in MV-CBCT multi-layer imager design.

Megavoltage (MV) cone-beam computed tomography (CBCT) using an electronic portal imaging (EPID) offers advantageous features, including 3D mapping, treatment beam registration, high-z artifact suppression, and direct radiation dose calculation. Adoption has been slowed by image quality limitations and concerns about imaging dose. Developments in imager design, including pixelated scintillators, structured phosphors, inexpensive scintillation materials, and multi-layer imager (MLI) architecture have been explored to improve EPID image quality and reduce imaging dose. The present study employs a hybrid Monte Carlo and linear systems model to determine the effect of detector design elements, such as multi-layer architecture and scintillation materials. We follow metrics of image quality including modulation transfer function (MTF) and noise power spectrum (NPS) from projection images to 3D reconstructions to in-plane slices and apply a task based figure-of-merit, the ideal observer signal-to-noise ratio (d') to determine the effect of detector design on object detectability. Generally, detectability was limited by detector noise performance. Deploying an MLI imager with a single scintillation material for all layers yields improvement in noise performance and d' linear with the number of layers. In general, improving x-ray absorption using thicker scintillators results in improved DQE(0). However, if light yield is low, performance will be affected by electronic noise at relatively high doses, resulting in rapid image quality degradation. Maximizing image quality in a heterogenous MLI detector (i.e. multiple different scintillation materials) is most affected by limiting total noise. However, while a second-order effect, maximizing total spatial resolution of the MLI detector is a balance between the intensity contribution of each layer against its individual MTF. So, while a thinner scintillator may yield a maximal individual-layer MTF, its quantum efficiency will be relatively low in comparison to a thicker scintillator and thus, intensity contribution may be insufficient to noticeably improve the total detector MTF.

Multi-layer imager design for mega-voltage spectral imaging.

The architecture of multi-layer imagers (MLIs) can be exploited to provide megavoltage spectral imaging (MVSPI) for specific imaging tasks. In the current work, we investigated bone suppression and gold fiducial contrast enhancement as two clinical tasks which could be improved with spectral imaging. A method based on analytical calculations that enables rapid investigation of MLI component materials and thicknesses was developed and validated against Monte Carlo computations. The figure of merit for task-specific imaging performance was the contrast-to-noise ratio (CNR) of the gold fiducial when the CNR of bone was equal to zero after a weighted subtraction of the signals obtained from each MLI layer. Results demonstrated a sharp increase in the CNR of gold when the build-up component or scintillation materials and thicknesses were modified. The potential for low-cost, prompt implementation of specific modifications (e.g. composition of the build-up component) could accelerate clinical translation of MVSPI.

Spectral imaging using clinical megavoltage beams and a novel multi-layer imager.

We assess the feasibility of clinical megavoltage (MV) spectral imaging for material and bone separation with a novel multi-layer imager (MLI) prototype. The MLI provides higher detective quantum efficiency and lower noise than conventional electronic portal imagers. Simulated experiments were performed using a validated Monte Carlo model of the MLI to estimate energy absorption and energy separation between the MLI components. Material separation was evaluated experimentally using solid water and aluminum (Al), copper (Cu) and gold (Au) for 2.5 MV, 6 MV and 6 MV flattening filter free (FFF) clinical photon beams. An anthropomorphic phantom with implanted gold fiducials was utilized to further demonstrate bone/gold separation. Weighted subtraction imaging was employed for material and bone separation. The weighting factor (w) was iteratively estimated, with the optimal w value determined by minimization of the relative signal difference ([Formula: see text]) and signal-difference-to-noise ratio (SDNR) between material (or bone) and the background. Energy separation between layers of the MLI was mainly the result of beam hardening between components with an average energy separation between 34 and 47 keV depending on the x-ray beam energy. The minimum average energy of the detected spectrum in the phosphor layer was 123 keV in the top layer of the MLI with the 2.5 MV beam. The w values that minimized [Formula: see text] and SDNR for Al, Cu and Au were 0.89, 0.76 and 0.64 for 2.5 MV; for 6 MV FFF, w was 0.98, 0.93 and 0.77 respectively. Bone suppression in the anthropomorphic phantom resulted in improved visibility of the gold fiducials with the 2.5 MV beam. Optimization of the MLI design is required to achieve optimal separation at clinical MV beam energies.